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ABSTRACT 

Short-period binaries represent extreme cases in the generation of stellar coronae via a ro- 
tational dynamo. Such stars are important for probing the origin and nature of coronae in the 
regimes of rapid rotation and activity saturation. VW Cep (P = 0.28 d) is a relatively bright, 
partially eclipsing, and very active object. Light curves made from Chandra/HETGS data show 
flaring and rotational modulation, but no eclipses. Velocity modulation of emission lines indicates 
that one component dominates the X-ray emission. The emission measure is highly structured, 
having three peaks. Helium-like triplet lines give electron densities of about 3 — 18 x 10 10 cm~ 3 . 
We conclude that the corona is predominantly on the polar regions of the primary star and 
compact. 

Subject headings: Stars: coronae - Stars: X-rays - Stars: Individual, VW Cep 



1. Introduction 

VW Cep (HD 197433), one of the X-ray- 
brightest of contact binaries, is a W-type W UMa 
system — one in which the more massive and 
larger star has lower mean surface brightness such 
that the deeper photometric eclipse occurs during 
the occultation of the smaller star. VW Cep has 
an 0.28 d (24 ks) orbital period, is partially eclips- 
ing (i = 63°), and has component spectral types of 
K0 V and G5 V (Hill 1989; Hendry & Mochnacki 
2000), or, according to Khajavi, Edalati & Jassur 
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(2002), F5 and GO. 

Among the coronally active binaries, activity is 
strongly correlated with the rotation rate, and in 
particular, the Rossby number, which is a mea- 
sure of the relative importance of Coriolis forces 
in the convective layer, which is in turn related 
to the magnetic dynamo strength (Pallavicini 
1989). At periods below one day, activity satu- 
rates (Vilhu & Rucinski 1983; Cruddace & Dupree 
1984). Since the activity level, as defined by 
L x / Xboi, actually decreases with increasing rota- 
tion rate, this trend has been referred to as "super- 
saturation" (Prosser et al. 1996; Randich 1998; 
Jardine & Unruh 1999; James et al. 2000; St§pieri, 



Schmitt & Voges 2001). The saturation level of 
log(L x /Lboi) ~ —3 extends down to periods of 0.4 
days. For contact binaries with periods between 
0.2 to 0.4 days, the median L x jL\ 30 \ is lower by 
a factor of four. Buzasi (1997) showed that in 
rapidly rotating low-mass stars, magnetic loops 
would be swept to the poles. Super-saturation 
may occur because loops are large and are un- 
stable to the Coriolis forces as they exceed the 
co-rotation radius: extended loops get swept to 
the poles (Jardine & Unruh 1999). Or loops could 
be compact (relative to the stellar radius), and 
the dynamics of surface flows clear equatorial re- 
gions (Stepieri, Schmitt & Voges 2001). The two 
scenarios are similar in that activity is predomi- 
nantly polar, but they differ in an important re- 
spect: X-ray sources are either large volume and 
rarefied or low volume and dense. Assuming cor- 
relation between photosphcric spots and coronal 
emission, optical light curve modeling is consis- 
tent with either scenario: Doppler image maps 
of VW Cep (Hendry & Mochnacki 2000) showed 
large polar spots. X-ray light curve modeling has 
been difficult because of the high probability of 
confusion by flaring, lack of phase redundancy, 
and the unavailability of high-resolution spectral 
diagnostics. 

Previous high-energy observations have charac- 
terized the coronae of VW Cep. Choi & Dotani 
(1998) analyzed ASCA spectra, and found a flux 
of about 1 x 10" 11 ergs cm~ 2 s -1 , and two com- 
ponent model temperatures of 7 and 22 x 10 6 K 
(logT ~ 6.8 and 7.3) with about equal emission 
measures. They obtained significantly reduced 
abundances of Fe, Si, Mg, and O, but a Solar value 
for Ne. A flare also occurred during the observa- 
tion, with a factor of three increase in the count 
rate. They used the flare emission measure and 
loop-scaling models to estimate a density of about 
5 x 10 10 cm~ 3 . The Ginga observations (Tsuru 
ct al. 1992) showed a thermal plasma tempera- 
ture in excess of 10 8 K, a flux similar to that de- 
termined by Choi & Dotani (1998), no rotational 
modulation, and Fe K flux lower than expected. 

The value of L x /Lbo\ ~ —3.6 places VW Cep 
in the super-saturated regime (using the bolomet- 
ric value implied by St§piefi, Schmitt & Voges 
(2001)). 



2. Observations and Reduction of High 
Resolution X-Ray Spectra 

We observed VW Cep for 116 ks with the Chan- 
dra High Energy Transmission Grating Spectrom- 
eter (HETGS) on August 29-30, 2003 (observation 
identifier 3766). The instrument has a resolving 
power (E/AE) of up to 1000 and wavelength cov- 
erage from about 1.5 A to 26 A in two independent 
channels, the High Energy Grating (HEG), and 
Medium Energy Grating (MEG) . For more details 
on the HETGS, see Canizares et al. (2005). 

VW Cep has also been observed in X-rays with 
the Chandra Low Energy Transmission Grating 
Spectrometer (Hoogerwerf, Brickhouse & Dupree 
2003), and with the XMM-Newton X-ray obser- 
vatory (Gondoin 2004b). The only other contact 
binary observed with HETGS is 44 Boo (Brick- 
house, Dupree & Young 2001). Observations with 
the HETGS provide the improved spectral resolu- 
tion and sensitivity required to better determine 
line fluxes, wavelengths, and profiles. Ultimately, 
we wish to determine which coronal characteris- 
tics are truly dependent upon fundamental stellar 
parameters. Combination of VW Cep data with 
that from LETGS and XMM-Newton, as well as 
with results for other short period systems (e.g., 
44 Boo, ER Vul) will help us to understand the 
basic emission mechanisms. 

Data were calibrated with standard CIAO 
pipeline and response tools 2 (CIA0 3.2 Thursday, 
December 2, 2004 / ASCDS version number). 
The updated calibration database geometry file 3 
provided a significant improvement in wavelength 
scales which resulted in good agreement between 
HEG and MEG as well as between positive and 
negative diffraction orders. This was important 
for line centroid analysis. 

Line, emission measure, and temporal analy- 
ses were done in ISIS (Houck & Denicola 2000; 
Houck 2002) and with custom code written in the 
S-Lang scripting language 4 using ISIS as a de- 
velopment platform. Figure 1 shows the 
counts spectrum obtained from the HETGS full 
exposure. The spectrum is qualitatively typical of 
coronal sources: a variety of emission lines from 

2 http://cxc. harvard.edu/ciao 

3 telD1999-07-23geomN0005 .fits 

4 S-Lang is available from http://www.s-lang.org/ 
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Fig. 1 — HETGS spectrum of VW Cep, 116 ks exposure. MEG (16433 counts) is the black line, and HEG 
(5269 counts) the gray. Insets show detail of the HEG spectrum (left) and MEG O VII triplet and N VII 
region (right). 



highly ionized elements formed over a broad tem- 
perature region, from O VII, N VII, Ne IX, and 
Fe XVII (logT ~ 6.3-6.7), up to high temperature 
species like S XV, S XVI, Ca XIX, and Fe XXV 
(logT ~ 7.2-7.8). It is apparent that iron has 
a fairly low abundance relative to neon, given the 
relative weakness of the 15A and 17A Fe XVII lines 
relative to the Ne IX 13A lines. The observed flux 
in the 2-25 A is 8.4 x 10~ 12 ergs cm~ 2 s -1 , and the 
luminosity (for a distance of 27.65 pc (Perryman 
et al. 1997) is 7.7 x 10 29 ergss" 1 . 

3. Analysis 

3.1. Light and Phase Curves 

The exposure lasted for five revolutions with- 
out interruption. Such phase redundancy is im- 
portant to discriminate intrinsic variability from 
that caused by rotational or eclipse modulation 
(for example, see the UV and X-ray monitoring 
studies of AR Lac, Neff et al. (1989); Huenemo- 
erder et al. (2003b)). Figure 2 shows light curves 5 
of diffracted photons (HEG and MEG, orders —3 
to 3, excluding zeroth) covering the entire HETGS 
spectrum as well as for two wide bands covering 

5 Light and phase curves were made with custom software 
(the aglc S-Lang package) available on the Chandra 
Contributed Software site: 

http: //exc .harvard . edu/cont-sof t/sof tware/aglc .1.2. 3 .html 
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Fig. 2. — Count rates of VW Cep in 2 ks bins. In 
the upper panel are light curves extracted from all 
diffracted photons (1.5-26 A; upper heavy curve), 
a "soft" band (12-26A; lower, thinner curve), and 
a "hard" band (1.5-8.3A; lower, thicker curve 
(blue)). The lower panel shows a hardness ratio 
(solid histogram) , and the median of the hardness 
for the first 70 ks (light dotted line). We defined 
the non-flare state as the first 80 ks. At the top 
we show an axis in which the integer part gives 
the number of rotations and the fractional part is 
the phase. 
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short ("hard", 1.5-8. 3 A) and long ("soft", 12- 
26 A) wavelength regions. The hardness counts 
ratio (defined as (hard — soft) /(hard + soft)) 
shows that the large increase after 80 ks is prob- 
ably a flare: proportionally more flux was emitted 
at shorter wavelengths, which are very sensitive 
to high temperatures through thermal continuum 
emission. Conversely, the bump in count rate be- 
tween 20-30 ks does not show in hardness, and 
likely is due to rotational modulation. 

We assume that flares are hot and will thus 
show a change in hardness, and that changes in 
volume alone will primarily show a change in count 
rate. More complicated scenarios are possible, 
such as extended loops with temperature gradi- 
ents rotating in and out of view and causing both 
hardness and count-rate variability on the time 
scale of a rotation. However, we know that flares 

- defined as hot, impulsive events with a grad- 
ual decay, in analogy to flares resolved on the Sun 

- are frequent on most active stars, and we will 
adopt the assumption as our working hypothesis. 
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Fig. 3. — The dark, solid histogram is the 
phase-folded X-ray count-rate curve of VW Cep, 
excluding flare times. Bins are A<f> = 0.025, 
or 600 s, with a cumulative exposure per bin 
of 1800 to 2400 s, depending on the phase. 
Errorbars on the histogram are the la statistical 
uncertainties. The gray points are a scaled 
visual intensity light curve to show the very 
different character of the optical and X-ray 
light spatial distributions (Visual photometric 
data are with permission from T. Pribulla at 
http : //www . astro . sk/^pribulla/lc/vwcep . dat .) 

Rotational modulation is apparent in the phase- 
folded curve of the non-flare times, which cov- 



ers about three orbital periods. Figure 3 shows 
this curve, using the ephemeris of Pribulla et al. 
(2002). No eclipses are obvious but there is grad- 
ual modulation on the orbital period with an am- 
plitude of 20%. The difference between the max- 
imum and minimum, given the la statistical er- 
rors shown, is about 5a. The visual light curve 
(Pribulla, Parimucha & Vanko 2000) has a sim- 
ilar amplitude but is very different qualitatively, 
having strong minima at phases 0.0 and 0.5, and 
continuous variability in between (a trademark of 
W UMa systems). The optical phase modulation 
is primarily dependent upon the system geometry 
and inclination. 

Without additional information, the X-ray light 
curve is not sufficient for localizing emission to one 
star or another. All we can say is that there is 
some asymmetric distribution, and possible occul- 
tation of longitudinally extended structures. It 
does not appear that there is significant emission 
from the smaller star, since there is no secondary 
eclipse at phase 0.0. The X-ray light has a distri- 
bution very different from the optical light. 

3.2. Velocity Modulation 

In the high-resolution spectrum, additional in- 
formation is available in the line positions which 
can help localize the emission by means of a de- 
termination of the mean radial velocity of the 
emitting plasma. In any single feature, even the 
strongest, Ne X (12.1 A), the line position is 
poorly constrained in phase bins small enough to 
sample the orbital velocity. Huenemoerder & Hu- 
nacek (2005) 6 showed that the Ne X mean Doppler 
velocity followed the primary's orbital radial ve- 
locity except for a sharp rise and fall across the 
secondary eclipse. This "flip-flop" is known as 
the Rossitcr effect (Rossiter 1924) if it is due to 
resolution of the rotational velocity profile of the 
star through occultation of velocity ranges during 
a transit. Given the low signal level available in a 
single line, this interpretation was not firm. 

We have improved upon that measurement by 
accumulating signal from many lines and con- 
structing a composite line profile. This mixes res- 
olutions and line shape since long- wavelength lines 
have higher velocity resolution than shorter wave- 
lengths, but it does have the advantage of increas- 

Preprint: http://arxiv.org/abs/astro-ph/0409258 
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ing signal greatly. Since HEG and MEG have dif- 
ferent resolutions and coverage, we kept their com- 
posites separate. To avoid blending, we chose only 
fairly isolated features, and these are flagged in the 
"Use" column of Table 1 with an "H" or an "M" . 

We accumulated spectra into 24 phase bins for 
each of the HEG and MEG positive and negative 
first orders (96 distinct spectra). We then trans- 
formed about 15 lines in each spectrum from wave- 
length to a common velocity scale and summed 
them. Finally, we combined plus and minus orders 
over a range of phase bins (to further improve the 
signal) and fit Gaussians to the composite profiles' 
cores (defined to be where the counts were greater 
than the maximum divided by 2e) separately for 
each of HEG and MEG. This method is similar to 
that used by Hoogerwerf, Brickhouse & Mauche 
(2004). We found that using five phase bins (a 
running average over A<f> = 0.2) was adequate to 
provide enough phase resolution and signal with- 
out losing too much sensitivity to velocity changes. 
To quantify the significance of the fitted velocity 
centroid, we computed the 90% confidence inter- 
vals (1.6c) for the velocity. Figure 4 shows 
the results. The composite line profiles make 
the background image, with velocity on the y-axis, 
and phase on the x-axis, and the darker shading is 
for higher intensity as indicated by the color-bar 
to the right, mapping color to counts. The yel- 
low error bars are the 90% confidence limits of the 
composite profiles centroid averaged over 5 phase 
bins (Acj) = 0.2). The red curves show the stellar 
center of mass velocities, also averaged over the 
same phase range as the data. 

The composite velocity centroid very closely fol- 
lows the radial velocity of the primary (more mas- 
sive) star. The preliminary sharp velocity tran- 
sition at <fi = 0.5 seen in neon (Huencmoerder & 
Hunacck 2005) is not apparent in the MEG curve, 
but the background HEG image does have a hint 
of a transition sharper than that due to the stellar 
radial orbital velocities. 

There are significant systematic deviations from 
the orbital velocities apparent near cf> = 0.8. 
Given the statistical uncertainties shown are 90% 
(1.6ct), the deviations from the expected velocity 
are about 2.5<7 for each grating. 

Note that the instrumental resolution is approx- 
imately 500 km s^ 1 FWHM at 12A for the MEG, 
and about half that for the HEG. We are able 



to determine centroids to much higher precision 
(~ 25kms" 1 for MEG; ~ 45kms" 1 for HEG) be- 
cause have accumulated signal over multiple pro- 
files, and because each profile is over-sampled by 
the spacecraft dither which randomizes line place- 
ment with respect to pixel boundaries. Hoogerw- 
erf, Brickhouse & Mauche (2004) achieved preci- 
sion of about ~ 15 kins -1 from MEG spectra, and 
also quote a statistical verification of the precision 
of a Gaussian fitted centroid. Also for comparison, 
Ishibashi et al. (2006) 7 using different techniques 
on HETGS spectra of Capella, a low-orbital ve- 
locity sytem, obtained 3a velocity uncertainties of 
about 20kms _1 which is about 11 kms -1 for 90% 
confidence. It is not surprising that our precision 
is somewhat less given the lower signal per phase 
bin and given more velocity smearing in phase for 
this short period system. 

The implications of the velocity centroid varia- 
tions with regard to geometric distribution of the 
emitting plasma will be discussed in Section 4.1 in 
conjunction with light curve variability, emission 
measure, and density determinations. 

3.3. Line Fluxes 

Line fluxes, which are required for differen- 
tial emission measure (DEM) analysis, were mea- 
sured by fitting a the sum of a continuum model 
and a number of Gaussians to narrow regions of 
the counts spectra. The parameters of the fits were 
the Gaussian centroids and fluxes. The contin- 
uum and Gaussian widths were determined a pri- 
ori. For the continuum, we used a plasma model 
determined iteratively from the emission measure 
solution (see Section 3.4). For the first iteration, 
we fit a two-temperature plasma model to rela- 
tively line-free regions, as determined by the model 
and instrument response. Since lines were measur- 
ably broader than the instrumental resolution, we 
determined the amount of excess line width re- 
quired at one wavelength for a strong, relatively 
isolated line (Mg XII 8.4A), then scaled the width 
with wavelength, under the assumption that the 
broadening is due to cither rotational or orbital 
effects. However, we simply treated the broad- 
ening as a Gaussian, which was adequate for the 
purposes of obtaining good fits to line fluxes and 
positions. The excess required was equivalent to 

7 Preprint: http://arxiv.org/ abs/astro-ph/0605383 
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about 200 km s 1 of turbulent broadening. 

All fits were done by convolving the model with 
the instrument response, using effective areas and 
line-responses made with CIAO tools (mkgarf and 
mkgrmf ; see Section 2). Positive and negative or- 
ders were combined dynamically (that is, in mem- 
ory, with no new counts or response files created 
on disk), and the HEG and MEG spectra were fit 
jointly in regions were there were sufficient counts 
in each. 

Line identification and blending were deter- 
mined iteratively with the DEM solution. Given 
the plasma model, blended components were 
added and removed from the fits according to 
whether the model showed that resolved blends 
were important or negligible. 

The measurement methods used here are sim- 
ilar to those described by Huenemoerder et al. 
(2003b). All the line fitting was done with scripts 
programmed in ISIS version f .2.8. 

The resulting counts spectrum and convolved 
model are shown in detail in Figure 5. 
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Fig. 4. — Composite line profile for the MEG spec- 
trum (top) and HEG (bottom). Yellow bars are 
the 90% confidence intervals of the centroid of the 
composite profile. The red sinusoidal curves are 
the center-of-mass velocities of the binary compo- 
nents. The background is an intensity map of the 
composite profiles' counts. The bar to the right 
of each image gives the counts-to-color mapping 
for the image. The systematic deviations from the 
primary's (more massive star) velocity are clear 
in the HEG for phases 0.5-1.0. The centroids are 
correlated since they were measured in phase bins 
of width 0.2. The red curves were also averaged 
over the same phase intervals as the data. 
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Fig. 5. — Detailed spectrum and synthetic spectrum from the DEM and abundance fit. The inner two 
panels of each box show the MEG (upper) and HEG spectra. The counts are shown in dark color (blue), 
the folded model in lighter color (red). The overlap is filled with black. The very top and bottom panels of 
each box are counts residuals for each spectrum. Lines from Table 1 are marked. Spectra were binned by 
two and Gaussian smoothed, except for the 1.7-4.7 A region where both HEG and MEG were grouped by 
4 bins, the 4.6-7.6 A and > 13 A regions where HEG was grouped by 4 bins. 
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All the line measurements are given in Table 1. 
We list more features than were used in analy- 
sis. Some were rejected because they had large 
wavelength residuals and so are misidentified or 
are strongly blended. Some were rejected because 
they had large flux residuals and so are misiden- 
tified or have poor atomic data. For example, the 
Fe xvn 1 5 A line is nearly always under-predicted, 
which may be due to deficiencies in the atomic 
data (Laming et al. 2000; Doron & Bchar 2002; 
Gu 2002). Some lines were fit simply to provide a 
good determination of the flux and wavelength of 
a nearby "interesting" feature. The "Use" column 
of Table 1 indicates which lines were used in emis- 
sion measure reconstruction, or for composite line 
profiles. The density sensitive He-like triplet inter- 
combination and forbidden lines were not used in 
DEM reconstruction. 

3.4. Emission Measure 

The differential emission measure is a one di- 
mensional characterization of a plasma, and can 
be defined as N% dV/dT, in which N e is the elec- 
tron density, V is the emitting volume, and T the 
temperature. The DEM is an important quantity 
because it represents the radiative loss portion of 
the underlying heating mechanism. An emission 
measure can be derived from measurements of line 
fluxes and some assumptions about the homogene- 
ity and ionization balance of the emitting plasma. 
Derivation of the emission measure relies on de- 
tailed knowledge of fundamental atomic parame- 
ters. Even given accurate emissivities, the contri- 
bution functions versus temperature are broad so 
the emission integral cannot be formally inverted. 
Hence there are many methods to regularize the 
solution to obtain the emission measure and ele- 
mental abundances. 

Huenemoerder & Hunacek (2005) used a sim- 
ple method described by Pottasch (1963), in which 
the DEM is approximated by a ratio of line lumi- 
nosity to average line emissivity at the tempera- 
ture of the maximum emissivity. Here we improve 
upon that by simultaneously fitting the DEM and 
abundances using a method similar to that de- 
scribed by Huenemoerder et al. (2003b), who also 
discussed some of the caveats of emission measure 
reconstruction and gave relevant citations. Briefly, 
the relation between emissivity and flux is an ill- 
posed problem. To obtain a unique solution, under 



assumptions of the model, a regularization term 
of some form must be included. Here we have 
replaced explicit smoothing of the DEM in our 
prior work with a regularization term, so that the 
functional form of the statistic we minimize is 

X 2 = £4 [Li A m £ t S t e lt D t ] 2 + qP(D)(l) 

in which I is a spectral feature index and t is the 
temperature index. The measured quantities are 
the line luminosities, L;, with uncertainties 07. 
The a priori given information are the emissivities, 
eu as defined by Raymond & Brickhouse (1996), S t 
is the logarithmic bin size, and the source distance 
(which is implicit in L). The minimization pro- 
vides a solution for the differential emission mea- 
sure, D, and abundances of elements Z, Az- To 
naturally constrain the DEM to be positive, we 
actually fit \nDt. P(D) is a regularization term 
(or "penalty function" ) which is only a function of 
the model, and q is a scale factor which specifies 
the relative importance of the regularization. For 
P, we used the sum-squared second derivative of 
the DEM with respect to logT, and a multiplier 
to make the two terms of the statistic comparable. 
This form imposes a minimum smoothness on the 
solution; the DEM cannot have large changes in 
curvature. We did not impose any regularization 
on the abundances, since we have no intuitive bias 
on their functional form. 

The DEM fitting was done with custom soft- 
ware in ISIS. While ISIS has no built-in DEM 
reconstruction, it does have sophisticated plasma 
database access and evaluation functions permit- 
ting efficient lookup and computation of emissivi- 
ties. It also provides for user-defined fit functions 
and statistics, which made it straightforward to 
develop a DEM model within its fitting and mod- 
eling infrastructure. 

As an initial guess, we assumed cosmic abun- 
dances and a boxcar DEM; the amplitude of the 
central temperature region was chosen to approx- 
imately produce the observed line counts, and the 
ends were set to a very low value. The reason for 
leaving the DEM at the temperature extrema low 
is that the DEM is not well constrained phys- 
ically there by lines, but the hot portion does 
strongly affect the continuum. Hence, we pre- 
ferred a solution in which emission weights are in- 
creased only if required by the fit to lines alone. 
The normalizations outside our regime of temper- 
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ature sensitivity, roughly logT = 6.4 (O Vll) to 7.8 
(Fe xxv; the tail of Si xvi), were artificially con- 
strained to have negligible emissivity. The line- 
to-continuum ratio was fit post facto by scaling the 
normalization of the DEM and abundances. 

The model DEM was obtained iteratively. For 
each trial DEM , we generated a new continuum 
model and re-fit the lines. After the second such 
iteration, we reviewed line residuals and excluded 
those which were very large (greater than 3c) , un- 
der the assumption that they were misidcntified or 
blended. The final fit included only the selected 
lines, and was post facto renormalized to repro- 
duce the observed line-to-continuum ratios in the 
8-11 A range. The line fluxes predicted by the 
DEM and abundance model are listed in Table 1, 
along with the residuals from the parametric fit, 
and their 8\ value (last three columns). The 
"Use" column flags lines used in the DEM recon- 
struction with an "E" . Due to the iterative nature 
of the model, the measured line fluxes are slightly 
dependent upon the DEM solution through the 
definition of the continuum. 

Figure 7 shows our final emission measure (inte- 
grated over logarithmic bins of width 0.1 dcx), us- 
ing the Astrophysical Plasma Emission Database 
(APED, Version 1.3.1) for line emissivities (Smith 
et al. 2001), the ionization balance of Mazzotta 
et al. (1998), and Solar abundances from Anders 
& Grevesse (1989). The envelope shown was de- 
termined by computing the deviation in solutions 
over a Monte-Carlo run of about 100 fits in which 
the line fluxes were varied according to their mea- 
sured statistical uncertainties, assuming a Gaus- 
sian distribution about their fitted values. This 
assumption would actually result in an enlarged 
uncertainty bounds, since the data are an esti- 
mate of the true mean, and the random pertur- 
bation would sometimes be toward the true mean, 
and sometimes away, perhaps further than the the 
measured uncertainty and truth would allow. 

Other sources of systematic error, the uncer- 
tainties in the atomic physics and in the instru- 
mental calibration, were not included, and proba- 
bly accounts for a similarly sized envelope. While 
systematic uncertainties in atomic data are diffi- 
cult to incorporate explicitly, Huenemocrder et al. 
(2003b) applied an approximate method, in which 
they set a global lower limit on the line flux un- 
certainty of 25% (regardless of the statistical un- 



certainty) and repeated the DEM reconstruction. 
They found a similar shaped distribution, but 
proportionally larger uncertainty. The important 
quantity for understanding coronal activity is the 
overall shape of the DEM, which seems reliable 
under the current assumptions and calibration un- 
certainties. 

There is sharp structure in the DEM with a 
large peak at log T — 6.9, a second peak at 6.6, and 
hot tail from about 7.2 to 7.8. This is qualitatively 
similar to the simple provisional DEM of Huen- 
emoerder & Hunacek (2005), but with sharper 
structure as expected from an iterative solution. 

The fitted coronal abundances, relative to So- 
lar (since we do not know the stellar photosphcric 
values), are shown in Figure 8. The uncertain- 
ties are the la range as determined by the Monte- 
Carlo iteration, except for Ni and Ca. The latter 
two were not included in the DEM solution since 
they were weak and were not fit with parametric 
profiles as were the strong lines (and so they do 
not appear in Table 1). Instead, we fit them post 
facto by using the DEM and abundance solution 
to define a plasma model, then adjusted the abun- 
dance of Ca by minimizing the residuals between 
the binned counts and synthetic spectrum in the 
Ca XIX He-like triplet near 3.2 A as a function of 
Ca abundance. There is a blend of Ar XVII here, 
but given the DEM , it is expected to be about an 
order of magnitude weaker than Ca. For Ni, we 
similarly fit the Ni XIX lines at 14.043 and 14.077 
A, and Ni XIX 12.435 A. In each region, blended 
Fe lines had model fluxes 5-10 times weaker than 
the nickel lines. The Ca abundance uncertainty is 
large because there are few counts and a relatively 
strong continuum. The Ni has smaller uncertainty 
because there are more counts and weaker contin- 
uum. 

3.5. Density 

The helium-like triplets are well known as den- 
sity sensitive diagnostics, particularly the ratio 
of the forbidden (/) to inter-combination (i) line 
fluxes, commonly known as the i?-ratio, R = f/i 
(Gabriel & Jordan 1969; Porquet & Dubau 2000; 
Porquct et al. 2001; Ness et al. 2001). A mildly 
temperature sensitive diagnostic is the G-ratio, 
which is the sum, f + i, divided by the resonance 
line (r) flux. The critical density increases with 
atomic number. In the HETGS range, the ions of 
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Fig. 7. — In this reconstructed emission measure 
distribution, the width of the filled region is the 
1-standard-deviation range as determined from a 
Monte-Carlo run in which the observed line fluxes 
were randomized according to their statistical un- 
certainty. The fit is physically unconstrained at 
the highest and lowest temperatures plotted, but 
was artificially constrained to have a negligible 
value at the extrema. The integrated emission 
measure is about 4.6 x 10 52 cm~ 3 . 




Fig. 8. — The reconstructed abundances relative 
to Solar are plotted against their first ionization 
potential (FIP). Error bars are the one standard 
deviation range from the Monte-Carlo emission 
measure and abundance reconstruction, or in the 
cases of Ni and Ca, from a post facto fit of abun- 
dances given the DEM. 



interest for coronal diagnostics are O VII, Ne IX, 
and Mg XL 

The O VII triplet ratios are relatively straight- 
forward to determine since the lines are well sepa- 
rated, unblended, and the continuum is low. Ne IX 
is a difficult case since it is blended with several 
line of Fe XIX and Fe XX, in addition to hav- 
ing a significant continuum. Ness et al. (2003a) 
performed a thorough analysis of this region in 
Capella at different resolutions, and showed that 
HEG resolution is necessary for accurate model- 
ing. The Mg XI region is of intermediate difficulty; 
it contains blends with the Ne X Lyman-like se- 
ries, which converges near the inter-combination 
line. Neon is often overabundant in stellar coro- 
nae, making the contribution possibly significant. 

We have determined R and G ratio confidence 
contours for O VII and Ne IX by directly fitting the 
line ratios to the spectra. Fits were done similarly 
to the lines (see Section 3.3), but with the addition 
of parameter functions to define fit-parameters for 
the ratios themselves. 

Figure 9 shows the result for oxygen. Here we 
have decided to use the LETGS data (observa- 
tion ID 2559), since it has better signal at O vn, 
and the contours are smaller than for the HETGS. 
Though at a different epoch, the solutions from 
HETGS and LETGS are equivalent. We have plot- 
ted the axes reversed, since the temperature and 
density increase inversely with the ratios. The 
contours show the one-, two-, and three-sigma con- 
fidence intervals, with the best fit marked with a 
"plus" sign. The underlying grid shows the the- 
oretical ratios' lines of constant temperature and 
density (density dependent tables are from Brick- 
house, private communication). We also show the 
counts spectrum and best fit model. The R ra- 
tio contours clearly show a density above the low- 
density limit; the best fit value is logN e = 10.5, 
with the lcr confidence interval ranges from 10.3 
to 10.8. The G-ratio spans a broad range but is 
not inconsistent with the reconstructed DEM. 

For the Ne IX region, we included nine resolved 
lines of Fe XVII, Fe XIX, and Fe XX. There is an 
unresolved blend of Fe XIX and Fe XX with the 
inter-combination line. From fits to the neighbor- 
ing lines, the DEM model and APED, we can 
estimate the contribution of iron to the Ne i-line 
of about 15%. This was included as a frozen 
component in the ratio fit. We show the counts, 
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Fig. 9.— The LETGS VII spectral region 
(top, dark histogram), model (top, light his- 
togram), and the confidence contours (bottom, 
colored ovals) for the G and R ratios. The the- 
oretical isothermal, isodensity G — R ratio grid is 
plotted in gray, and labeled with the logarithmic 
density and temperature values. The best fit den- 
sity is about log iV e = 10.5. The HETGS data give 
a similar result, but with larger contour regions. 
(The vertical spike in the contours is a numerical 
artifact.) 



Ne IX, and Fe model in Figure 10 for both HEG 
and MEG, since resolution is critical. The fig- 
ure also shows the model continuum, and Table 1 
flags the resolved iron blends included. For the 
unresolved blends, we give the identifying informa- 
tion from the APED, which labels lines uniquely 
by the element, ion, upper (u), and lower (I) lev- 
els: Fe XIX A13.551 (it = 65, I = 1); Fe XX 
A13.535 (u = 109,107, I = 7), A13. 558, 13.565 
(u = 110,109, I = 8). 

In Figure 11, we show the resulting G — R con- 
tours and spectra with best-fit model. Here, it 
is not clear that the density is above the criti- 
cal value; the best fit is \ogN e = 11.25, and the 
one-sigma range is 10.0-11.4. The G-ratio con- 
tours extend to low temperatures, but the range 
is not inconsistent with the DEM. We 
also fit the Mg XI ratios. The ratio and uncer- 
tainty gives R = 2.7 ± 0.5. The theoretical low- 
density ratio is 2.9. Hence, we are well below the 
limit of density sensitivity (the critical density is 
about logiV e = 13). We included the high-rt Ne X 
Lyman-like series and iron blends in the triplet re- 
gion (Testa, Drake & Peres 2004), which we show 
in Figure 12. 

4. Discussion 

There is no question that coronal activity sat- 
urates at short periods, and that super-saturation 
occurs at even shorter periods. Why saturation 
occurs is still an open question among hypotheses 
of dynamo suppression by tidal interactions, coro- 
nal stripping via Coriolis forces, filling factors near 
unity, or other effects. With high resolution X- 
ray spectra, we cannot yet solve this problem, but 
we can provide additional details from individual 
case studies. We can make several deduction from 
the VW Cep light curves, line velocities, emission 
measures, and density diagnostics. We will also 
consider the similar super-saturated W UMa sys- 
tem, 44 Boo. 8 

4.1. Coronal Geometry 

From visual inspection of the VW Cep X-ray 
light curve, we see that X-rays vary in brightness 
quasi-sinusoidally at the period of the binary with 



'44 Boo is also officially designated "i Boo" 
mistakenly called as "44i Boo." 
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Fig. 10.— Ne IX region for MEG (top) and HEG 
(bottom), showing the counts (dark histogram), 
the Ne IX components (red-filled Gaussians), the 
iron components (green-filled Gaussians) , and the 
continuum model (gray- filled region) . The relative 
intensity of the iron line blended with the Ne IX 
inter-combination line was determined from neigh- 
boring iron lines. Note that the counts scale is 
logarithmic. 



Fig. 11. — Ne IX spectral region (top, dark his- 
togram), model (top, light histogram), and the 
confidence contours (bottom, colored ovals) for 
the G and R ratios. The theoretical isothermal, 
isodensity G — R ratio grid is plotted in gray, 
and labeled with the logarithmic density and tem- 
perature values. The best-fit density is about 
logiV e = 11.25. 
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Fig. 12. — Mg XI region for the MEG, showing 
the counts (dark histogram), And individual ionic 
components as smooth curves: Mg XI (red), Ne X 
and Fe (green) 



a modulation of about 20%. The shape is very 
different from the optical curve, which is deter- 
mined primarily by the geometric configuration of 
two stars with nearly equal surface brightnesses 
(tidally distorted spheroids which undergo par- 
tial eclipses) and secondarily by inhomogeneous 
brightness distributions (gravity darkening and 
starspots). Since we do not see obvious eclipses 
in X-rays (which could have a phase duration of 
about 0.2), we infer that the regions obscured 
probably do not contribute significant X-ray emis- 
sion within our detection sensitivity. Since the in- 
clination is 63°, and the primary has about twice 
the radius of the secondary, the regions covered 
are roughly the southern hemispheres of each star. 
The north-polar region of each star is always visi- 
ble, roughly within co-latitude of 30°, so any emis- 
sion from this region on either star will contribute 
a constant emission component (within projected 
area affects). The non- uniform emitting material 
can thus be roughly constrained to latitudes of 
— 60°, in one hemisphere only, with more emis- 
sion measure seen near phase (larger star in 
front). 

If we now consider the X-ray composite radial 
velocity curve, we can rule out significant emis- 
sion from the secondary (smaller) star, at the limit 
of about 30%. The ratio of photospheric surface 
areas is about 2:1. If we had 30% of the X-ray 
line flux originating from the secondary, it would 
have been apparent in the composite line profiles, 



but we do not see it (see Figure 4). We explored 
the limits of the secondary's signature in the com- 
posite profiles through simulations. We simulated 
spectra using the same model emission measure 
for each star, but doing a weighted sum at each 
phase, with each spectrum redshifted according to 
the phase. We used the same exposure times from 
each bin of the observed phase curve as well as ones 
100 times larger. We considered relative weights, 
primary to secondary of 0.7:0.3 (the ratio of ar- 
eas), 0.8:0.2, and 1:0. After the simulated data 
were produced, the composite lines were formed 
and measured in exactly the same way as for the 
observation. 

For the 0.7:0.3 weighting, modulation in the 
composite centroid is reduced to an undetectable 
size. For 0.8:0.2, modulation is reduced by about 
50%, to about the maximum deviation seen in the 
HEG near phase 0.8. For all emission from the 
primary, the curve follows the primary radial ve- 
locity, as expected. Thus, we find most emission 
at most phases can be attributed to the primary. 
At some phases, approximately 0.75-1.0, we may 
have about 20% from the secondary. Since there 
is no apparent eclipse of the secondary (at phase 
0.0, the primary eclipse), the emission from the 
secondary must be at high latitudes on the trail- 
ing hemisphere, but not polar or we would have 
a larger velocity centroid perturbation near phase 
0.25. 

We can next apply emission measure and den- 
sity information to estimate coronal extent, under 
some simplifying assumptions about the emitting 
plasma's geometry. If we assume that the 20% 
light curve intensity modulation is caused by a 
random distribution of identical coronal structures 
of constant cross-section, then we can have about 
25 structures, based purely on counting statis- 
tics. Using the scaling relation of Huenemoerder, 
Canizares & Schulz (2001), an integrated emission 
measure of 4.6 x 10 52 cm -3 and a density range 
of logiV e = 10.5 to 11.25, and assuming all emis- 
sion comes from the primary of radius 0.9 i?©, then 
the structure height is about 0.06 to 0.2 stellar 
radii. Such a corona would thus be fairly com- 
pact, which is self consistent; if it were very ex- 
tended, flux modulation could be much less. A 
single such structure could have a height from 0.2 
to 0.6 stellar radii and could probably be placed to 
create the same intensity modulation. We prefer 
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the multiple structure hypothesis, but this is by 
no means a unique requirement of the data. If we 
spread the emitting volume more uniformly over a 
hemisphere, then for the above density range the 
fractional height is 0.02 to 0.03. We probably have 
a compact, polar corona. 

Gondoin (2004b) reported an X-ray eclipse seen 
XMM-Newton light curves obtained 10 months 
prior to the CXO observations and concluded that 
both components were X-ray emitters. If the dip 
in XMM-Newton rate were really an eclipse, then 
this would indicate rather large and rapid changes 
in coronal structure. The XMM-Newton observa- 
tion, however, did not cover a single continuous or- 
bit, and the dip is offset from phase 0.0, so the in- 
terpretation is not definitive. Also contrary to our 
results, Gondoin (2004b) derived loops quite large 
in comparison to the stellar radii (20-80%), but ad- 
mittedly by methods which are not self-consistent. 

4.2. Temperature Distribution, Heating, 
Opacity 

The differential emission measure (Figure 7) is 
similar to other active stars, being highly struc- 
tured and spanning a broad range in temperature. 
Sanz-Forcada, Brickhouse & Dupree (2002, 2003) 
show a large collection of DEM distribtions with 
a variety of peaks, bumps, and tails. In our analy- 
sis, we included the flare times and this no doubt 
produces the bump above logT = 7.1, similarly to 
the results for II Peg (Huenemoerder, Canizares & 
Schulz 2001) and Proxima Cen (Giidel et al. 2004), 
which were modeled at different activity levels. 

The most prominent peak in the VW Cep 
DEM is well defined and narrow, with a slope ap- 
proximately proportional to T 4 . Such a feature 
has also been observed in other stars of differ- 
ent types (Brickhouse et al. 2000; Sanz-Forcada, 
Brickhouse & Dupree 2002). The DEM can give 
us some insights into the underlying physics. A 
possible interpretation of this recurrent feature is 
linked to the spatial location and temporal distri- 
bution of coronal heating. Testa, Peres & Reale 
(2005) have proposed a hydrodynamic model of 
loops undergoing pulsed heating at their foot- 
points, which is able to reproduce the presence 
of a peak and the steep rise of the DEM observed 
in these active stars. Whatever the mechanism, 
the heating in VW Cep does not look radically 
different from other coronal sources, at least as 



manifested in the DEM. 

The abundances do not show any strong trend 
with the first ionization potential (FIP). While Ar 
looks unusually abundant in Figure 8, the statis- 
tical uncertainty is large. Furthermore, the mea- 
surement is also subject to systematic uncertainty 
in placement of the continuum. The Ar abun- 
dance is probably not significantly different from 
Ne. The overall distribution looks very similar, for 
example, to that of the RS CVn binary, AR Lac 
(Huenemoerder et al. 2003a). The Ne:0 ratio 
which we derived from the iterative emission mea- 
sure and abundance reconstruction is 0.35 ± 0.03. 
This is comparable to the rather robust active 
star mean of 0.41 found by Drake & Testa (2005) 
from temperature-insensitive line ratios. From the 
abundance analysis, we again have no reason to 
think that the corona of a W UMa system, though 
super-saturated, differs from other coronally ac- 
tive stars. 

We examined line ratios which have been used 
for opacity diagnostics. Fe xvn has been problem- 
atic, because the ratios systematically differ from 
theoretical values (Laming et al. 2000; Doron & 
Behar 2002; Gu 2002). VW Cep is no exception; 
its Fe xvn 15.01 A line is over-predicted. We find 
no significant differences in ratios from typical op- 
tically thin values shown in Testa et al. (2004) or 
Ness et al. (2003b). 

4.3. Comparison to 44 Boo 

44 Boo is a contact binary system which is 
very similar to VW Cep, having similar period, 
mass ratio, and spectral types. Its fundamental 
characteristics are given by Hill, Fisher & Holm- 
gren (1989) and a velocity curve by Lu, Rucinski 
& Ogloza (2001). Brickhouse, Dupree & Young 
(2001) have analyzed HETGS spectra and con- 
cluded that most emission is at high latitudes, 
compact, that there is a very compact emitting 
region on the smaller star (the secondary) which 
causes a very brief dip in X-ray brightness when 
it momentarily rotates out of view, and a larger 
emitting region near the pole of the primary which 
gives rise to a quasi-sinusoidal brightness modula- 
tion, but which isn't eclipsed. Gondoin (2004a) 
presented XMM-Newton observation of 44 Boo 
which covered on orbital period. Based on the 
absence of eclipses, he concluded that the corona 
must be extended, though he presented other ar- 
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guments which led to compact loops. 

Using the same techniques applied to VW Cep, 
we have extracted the X-ray light and phase curves 
and a composite line velocity curve from the same 
HETGS observation used by Brickhouse, Dupree 
& Young (2001) (Observation ID 14). The light 
curve in Figure 13 shows three flares. The rise 
in the final 10 ks, identified as a flare by Brick- 
house, Dupree & Young (2001), shows no in- 
crease in hardness, so we accept it into the phased 
curve. Our phased X-ray light curve (Figure 14) 
is qualitatively different from Brickhouse, Dupree 
& Young (2001). 

Since W UMa stars undergo period variations, 
we were careful to reference the ephemeris to light- 
curves from the same epoch (Pribulla et al. 2001). 
We also obtained a B photometric light curve from 
the same epoch and scaled it to arbitrary relative 
intensity to display over the X-ray curve. We do 
find strong evidence for the primary eclipse in the 
broad dip in X-ray flux in good correspondence 
with the B intensity. The rest of the curve is qual- 
itatively similar to the VW Cep X-ray curve (Fig- 
ure 3), having a rather broad depression with no 
obvious eclipse. The gap beginning near phase is 
due missing data due to the exclusion of flares. We 
do not find a very narrow dip nor require a very 
small near-polar emitting region as hypothesized 
by Brickhouse, Dupree & Young (2001). These dif- 
ferences are probably primarily due to our forming 
the histogram in phase bins, weighting the counts 
by the exposure function (which ranges from zero 
to 2 ks), and in our flare discrimination using a 
hardness ratio. 

Figure 15 shows our composite line profiles 
and centroid, which is qualitatively similar to the 
velocity curve of Brickhouse, Dupree & Young 
(2001). The primary clearly dominates the 
composite line centroid. However, the distortions 
and wings also hint at some contribution from 
the secondary, particularly between phases 0.6- 
0.9 where the centroid is systematically shifted to 
more positive velocities. The secondary cannot 
contribute as much as half the light, however, or 
the composite profile centroid would be strongly 
biased toward zero. Given the presence of the pri- 
mary eclipse, the broadening in the profile wings, 
slight distortions in the HEG centroid, and simula- 
tions, we estimate that the secondary contributes 
no more than about 20% of the X-ray flux during 
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Fig. 13. — Count rates of 44 Boo in 200 s bins. In 
the upper panel are light curves extracted from all 
diffracted photons (1.5-26 A; upper heavy curve), 
a "soft" band (12-26A; lower, thinner curve), and 
a "hard" band (1.5-8.3A; lowest, thicker curve). 
Near the top is an axis giving the orbital phase 
plus a rotation count. The lower panel shows a 
hardness ratio (solid histogram), and the median 
of the hardness for the middle flat section (light 
dotted line). We defined the three high hardness 
regions as flares and excluded them from the phase 
curve (0-5, 12-18, and 38-46 ks). 



Fig. 14. — The dark, solid histogram is the phase- 
folded X-ray count-rate curve of 44 Boo, excluding 
flare times. Bins are A(j> = 0.025, or 462 s, with a 
cumulative exposure per bin of about 460 to 1400 
s, depending on the phase. The gap near phase 0.4 
is from missing data, due to flare-time filtering. 
The smooth curve is a scaled B light curve to 
compare optical and X-ray light spatial distribu- 
tions. Photometric B data are from T. Pribulla 
(http : //www . astro . sk/~pribulla/lc/ 44boo . jpg.) 
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Fig. 15. — Composite line profile for the MEG 
spectrum (top) and HEG (bottom) for 44 Boo. 
Yellow bars are the 90% confidence intervals of 
the centroid of the composite profile. The red si- 
nusoidal curves are the center-of-mass velocities of 
the binary components. The background is an in- 
tensity map of the composite profiles' counts. The 
bar to the right of each image gives the counts-to- 
color mapping for the image. The centroids are 
correlated since they were measured in phase bins 
of width 0.2. The red curves were averaged over 
the same intervals. 



some phases. 

Were it not for the primary eclipse, the 44 Boo 
situation would be very much like VW Cep: quasi- 
sinusoidal light curve with a minimum near phase 
0.5, 20% contribution from the secondary at some 
phases. 44 Boo is at a slightly higher inclination 
than VW Cep, 73° vs. 63°. This does give it 
a higher probability of eclipsing compact coronal 
emission concentrated at high latitudes. On the 
other hand, it is highly possible that the size and 
location of the coronal emission migrates, as do 
the photospheric spots (see, for example, Hendry 
& Mochnacki (2000)). At future epochs, it is pos- 
sible that VW Cep will show X-ray eclipses. 

5. Conclusions 

Super-saturation, at least in the cases of 
VW Cep and 44 Boo, appears to be manifested 
in small area of coverage by compact, near-polar 
structures. Given roughly a factor of two for coro- 
nal structures occurring in only one hemisphere, 
and a factor of two in expected emission per area 
(the secondary is weak relative to the primary), 
the X-ray flux can be significantly depressed below 
the saturated limit of L x /Lboi ~ 10 . Why this 
occurs is still unexplained, but we favor the gen- 
eral scenario of (Stejneh, Schmitt & Voges 2001) 
which forms polar, but compact, structures. 

Given the small area of corona on the sec- 
ondary, we predict that the secondary's coronal 
signature be highly changeable with phase over 
different epochs. This would show as changes in 
the depth of primary eclipse and phase of pertur- 
bation in the velocity profile. The HETGS is the 
only X-ray spectrograph currently capable of mak- 
ing these measurements. Obtaining high signal-to- 
noise per phase bin is of utmost importance to sur- 
face geometry reconstruction techniques. Further 
multi-orbit spectroscopy of these or other short 
period systems is certainly important to determi- 
nation of transient or common features. A factor 
of 10 increase in exposure would permit line pro- 
file diagnostics in a single line (Ne x) and greatly 
facilitate modeling. 

This research was supported by NASA grant 
G03-4005A and by NASA through the Smithso- 
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SV3-73016 for the Chandra X-Ray Center and Sci- 
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Table 1 — Continued 
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NOTE. — Values given in parentheses are the one standard deviation uncertainties on the preceding quantity. If 
the uncertainty has a value of " — ", then either the confidence did not converge, or the parameter was frozen. 
"Unused" features were fit in order to obtain a good fit to a region and determine values for nearby interesting 
lines. 

a Thc mnemonic is a convenience for uniquely naming each feature. It is comprised of the element and ion (in 
Arabic numerals) followed by a string indicating a wavelength and the wavelength (e.g., wl6.78), or a code for 
the hydrogen-like ("H") and helium-like "He" series, "L" for Lyman transition, one of "a", "b" , "g" , "d" , or "e" 
for series lines a, (3. 7, 8, e, and "r" , "i" , or "f" for resonance, intersystcm, or forbidden lines. 

b "Use" indicates whether the line was used in the emission measure reconstruction ("-" for no, "E" for yes); 
whether the feature was used in the composite line profile ("H" and "M" for HEG and MEG); or whether the 
line was used in the Nc IX triplet fit with the character, "t" . 

c Averagc logarithmic temperature [Kelvins] of formation, defined as the first moment of the emissivity distri- 
bution. 

d Theorctical wavelengths of identification (from APED), in A. If the line is a multiplet, we give the wavelength 
of the strongest component. 

c Mcasured wavelength, in A (uncertainty is in mA). 
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f Emitted source line flux is 10 -6 times the tabulated value in [phot cm -2 s -1 ]. 
s Model line flux is 10 — 6 times the tabulated value in [phot cm - 2 s~ 1 ] . 

Line flux residual, Sf — f a — ft- 
i Sx = Uo-ft)/<To- 

'The Ne IX intercombination line flux is blended with Fe XIX 13.551 A, and with Fe XX 13.535, 13.558, 13.565 
A lines. The tabulated flux includes these blends. About 85% of the flux is from neon, as determined by the 
neighboring lines. 

k The O VIII Lyman beta-like line is blended with Fe XVIII 16.004 A, whose flux is about 1.17 times that of 
Fe XVIII A 15.62, according to our DEM and APED model. Hence, we estimate that the actual 08HLb flux is 
45.9 (8.6). 
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